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ABSTRACT 

Context. Infrared (IR) spectra provide a prime tool to study the characteristics of polycyclic aromatic hydrocarbon (PAH) molecules 
in regions of star formation. Herbig Ae/Be stars are a class of young pre-main sequence stellar objects of intermediate mass. They are 
known to have varying amounts of natal cloud material still present in their direct vicinity. 

Aims. We characterise the IR emission bands, due to fluorescence by PAH molecules, in the spectra of Herbig Ae/Be stars and link 
observed variations to spatial aspects of the mid-IR emission. 

Methods. We analysed two PAH dominated spectra from a sample of 15 Herbig Ae/Be stars observed with the Spitzer Space 
Telescope. 

Results. We derived profiles of the major PAH bands by subtracting appropriate continua. The shape and the measured band char- 
acteristics show pronounced variations between the two Spitzer spectra investigated. Those variations parallel those found between 
three infrared space observatory (ISO) spectra of other, well-studied, Herbig Ae/Be stars. The derived profiles are compared to those 
from a broad sample of sources, including reflection nebulae, planetary nebulae, Hn regions, young stellar objects, evolved stars 
and galaxies. The Spitzer and ISO spectra exhibit characteristics commonly interpreted respectively as interstellar matter-like (ISM), 
non-ISM-like, or a combination of the two. 

Conclusions. We argue that the PAH emission detected from the sources exhibiting a combination of ISM-like and non-ISM-like 
characteristics indicates the presence of two dissimilar, spatially separated, PAH families. As the shape of the individual PAH band 
profiles reflects the composition of the PAH molecules involved, this demonstrates that PAHs in subsequent, evolutionary linked 
stages of star formation are different from those in the general ISM, implying active chemistry. None of the detected PAH emission 
can be associated with the (unresolved) disk and is thus associated with the circumstellar (natal) cloud. This implies that chemical 
changes may already occur in the (collapsing?) natal cloud and not necessarily in the disk. 
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1 . Introduction lution (0.5 - 10 Myr), the increased stellar winds from the cen- 

, 1 ttal object will have removed most of the surrounding molecular 

Herbig Ae/Be stars are a class of young stellar objects <|Herbig| doud materiaL leayi Qnl the KmnanU ive> acc retion disk 

[I960!) originally limited to pre-main sequence stars of interme- ^ ^oOALT M disk ~ 5.6 ■ KT 4 - 1.5 ■ 10 1 M ; HabartetalJ 

diate mass (2 - 8 M ) with spectral class A or B with associated 2004i Acke et al 2QQ ^ Waelkens & w^te7s][T^7h . This disk - 

nebulosity and showing emission lines. Over the years, however, tQ bfi ^ sifc of ^ formation . shows up m polarised 

the definition has been somewhat relaxed, e.g. the necessity for u ht and 

as a strong excess in infrared spectra. 

associated nebulosity (IMalfait et al.lll998l) . 

The formation and evo lution of isolated , low-mass stars, is With the 1995 launch of the Euro pean Space Agenc y's 

relatively well understood dShu et al.lll987l) . However, the gen- Infrared Space Observatory (ESA ISO, iKessler et~aT1.ll996l) . it 

esis of intermediate and high-mass stars remains more elusive, became possible for the first time to study the infrared spec- 

For intermediate mass, Herbig Ae/Be stars, it is thought that the trum of Herbig Ae/Be stars in detail. Those studies revealed a 

formation scheme is more or less a scaled up version of that for large variety in dust composition, both carbonaceous and sil- 



low- mass stars. The picture that has thus emerged is that stars icate in nature, and in dust properties (Ivan den Ancker et alj 



form in large molecular clouds and through successive stages 2000; lBouwman et al.ll2000t IWaelkens etaill 19961) . Specifically, 

of fragmentation, gravitational collapse, and disk accretion, ac- strong infrared (IR) features with satellite features and emission 

companied by bipolar outflows, the star reaches the main se- plateaus around 3.3, 6.2, 7.6, 7.8, 8.6, 11.2, and 12.7 micron, as- 

quence. During the later stages of the pre-main sequence evo- sociated with polycyclic aromatic hydrocarbons (PAHs), were 

apparent in many of these objects. PAHs, large molecules of 

Send offprint requests to: C. Boersma AT astro.rug.nl many fused aromatic rings, fluoresce in the infrared upon the ab- 
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sorption of a single visible or UV photon dLeger & P uget 1984; 
ICohen et alJll985t lAllamandola et al.lll989l) . Driven by the dis- 
persal of the disk and the formation of a planetary system, pro- 
toplanetary disks are continuously evolving. Geometry, gas, and 
dust content; gas and dust processing; and transport and mix- 
ing are important issues concerning these disks. The diversity of 
spectral features and their variations within and from source-to- 
source reflect these issues. However, the interrelationship of the 
variance in dust composition, the characteristics of the disk, and 
the global system is not well understood. 

Over the years, PAHs have proven to be an ideal probe 
for m any of these properties, e.g., disk geometry ( Meeu set al.l 
1200 lb . The presence of (large) PAH molecules influences many 
aspects in protoplanetary systems. These aspects include: (sur- 
face) chemistry, due to the large surface area of PAHs; heating 
and cooling through photo-electric ejection; infrared emission; 
gas-grain collisions; and the charge balance, which in its turn 
influences the equilibrium state of chemical reactions. 

Today, the superior sens itivity of NASA's Spitzer Space 
Telescope dWerner e t al. 2004) allows us to investigate proto- 
planetary disk systems in better detail than before. Spectra of 15 
nearby Herbig Ae/Be stars (Sect. GO), ob tained with the Infrared 
Spectrograph,] (IRS, iHouck et alTl2004l) . are investigated for 
emission due to PAHs (Sect. [3). We characterise the individual 
profiles for the two PAH-dominated spectra found in the sample 
(Sects. [3T1 & I3.2I I and, by using parallels between the morphol- 
ogy of the stellar systems of these sources and those of three 
well-studied systems, explore possible causes for the observed 
variations in band profiles (SecHJl. Similarities between the 6.2 
and "7.7" pm band profiles are treated (Sect. 14. It and the ex- 
tended nature of the PAH emission is investigated (Sect. |4.2t . 
The origin of the variations in peak position is discussed (Sect. 
14. 3t and from this, implications about the onset of the chemical 
modification of PAHs - disk and/or envelope - are inferred (Sect. 
14.4) . We finish by summarising our findings (Sect. [5J. 



2. Data 

The IRS instrument on boa rd Spitzer w as used to obtain the 
spectra of 36 Herbig Ae/Be (Herbig 1963) stars. W e selected the 
stars f rom a list of Herbig Ae/Be stars compiled by Malfait et al. 
(1998) and they complement the 26 Herbig Ae/Be stars ob- 
served in the Spitzer guaranteed time observations and legacy 
programs. 

For 15 genuine Herbig Ae/Be stars, we have access to the ap- 
propriate wavelength region for detecting emission due to PAHs. 
Those stars were observed using the short-high and short-low 
(SH, SL resp.) modules in high accuracy peak-up mode (<x ~ 
0.14"), providing spectra at R ~ 600 for the high-resolution 
module and R ~ 64 - 128 for the low -resolution module, cov- 
ering wavelengths from 9.9 - 19.6 and 5.2 - 14.5 pm for the SH 
and SL modules, respectively. 

The Spitzer Science Center (SSC) processed the data with 
the science data pipeline version S13.2.0. In order to make more 
accurate corrections for background emission than possible with 
the standard pipeline, 

we adopted the extraction techniques developed by the for- 
mation and evolution of protoplanetary systems (FEPS) and 
cores to disks (C2D) Spitzer legacy teams (Bouwman in prep., 



lLahuisll2007l and c2d-IRS delivery documentation^). Extraction 
of the one-dimensional spectra are done on the intermediate 
droopres and rsc products that are provided by the SSC. These 
data have most of the Spitzer-specific artifacts removed. For the 
rsc products this includes straylight and crosstalk corrections, 
but not flat-fielding. The absolute flux calibration is done using 
dedicated spectral response functions, derived from IRS spectra 
of a suite of calibrators using Cohen and MARCS stellar models 
provided by the SSC. 

The high-resolution spectra are extracted from the rsc prod- 
ucts using an optimal source profile extraction developed by 
the c2d team. The optimal extraction uses an analytical point 
spread function (PSF) defined from high S /N sky corrected cal- 
ibrator data. A simultaneous cross-disperion source profile and 
extended emission fit is made to the combined dither data. For 
the extended component, the cross-dispersion profiles of the flat- 
field images are used. The applied source profile is optimised 
by comparing the normalised (in cross-dispersion direction) and 
collapsed (in dispersion/wavelength direction) science source 
profile with the average profile of the calibrator (point) sources. 
Bad-pixels are identified and excluded from the profile fits. The 
optimal extraction allows us to correct the spectrum for local ex - 
tended emission (see, e.g., iGeers et alJ2006b l Lahu is et a 1. 2007), 
which is crucial in confused regions where the use of separate 
sky observations is limited. 

The low-resolution spectra are extracted from the droopres 
products using a 6.0 pixel fixed-width aperture. Background cor- 
rections are made by subtracting the two-dimensional spectral 
images from the associated nod positions. Besides stray light, 
this also corrects for anomalous dark current effects. Identified 
bad pixels are interpolated from the surrounding perimeter. The 
optimal source extraction is used to derive a possible resid- 
ual extended emission component in the extraction aperture of 
the low-resolution order 1 spectra (see Fig. [8] Due to the se- 
vere undersampling of the spectral images in the spatial di- 
rection, no source profile fits are made for order 2 and 3 and 
hence no estimate of the extended emission is derived. Flux 
offsets between orders (particular in the SL module) and mod- 
ules are apparent. They are probably related to pointing errors 
and are likely to have introduced low level fringing and small 
slope changes. For the SL module, we developed a procedure 
to estimate the pointing errors and correct for them. This pro- 
cedure minimises the total variance between the flux levels of 
the individual cycles, nod positions and orders, applying wave- 
length dependent correction factors to the flux for each relative 
offset. The procedure simultaneously corrects for fringes. For 
the high-resolut ion modules, fringes ar e removed using the irs- 
fringe package dLahuis & Boogertll2003l) . The SMART0 A mis- 
match between flux levels of the SL and SH modules is cor- 
rected by multiplying the SH flux levels with a constant, which 
is determined in the overlap region and no adjustment for any 
differences in tilt for the two sub-spectra is made. As the final 
reduction steps, outliers are removed and, in regions of spec- 
tral overlap, the orders are spliced to form a continuous profile. 
For rising spectra this splicing may introduce small artifacts in 
the overlap regions, most notably near 7.5 pm. This is probably 
related to weak order leakage effects introduced by the spectral 



1 The IRS was a collaborative venture between Cornell University 
and Ball Aerospace Corporation funded by NASA through the Jet 
Propulsion Laboratory and Ames Research Center. 



2 http://data.spitzer.caltech.edu/popular/c2d/ 

2006 1 20 1 _enhanced_v 1 /Documents/C2D_Spec tr_Expl_Supp.pdf 

3 SMART was developed by the IRS Team at Cornell University 
and is available th rough the Spitzer Science Center at Caltech. 
(Higdon et al. 2004) software package has been utilised to further re- 
duce and analyse the data. 
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response function, which has been determined using only declin- 
ing stellar spectra. In the case of extended sources, artifacts may 
arise due to small differences in aperture size between orders 
and modules. The relative errors are dominated by the intrinsic 
(Gaussian) noise i n each spectral d ata point, not by uncertainties 
in the calibration (Hine s et al . 2005). The intrinsic noise is char- 
acterised by the deviations found between subsequent nod posi- 
tions and cycles and are about 2 percent across an order. Between 
orders it is estimated to be about 5 percent and between modules 
about 8 percent. 

3. Analysis 

In order to investigate the emission due to PAHs we roughly 
classify the spectra of the 15 sources by their dominant spectral 
component. Table Q]presents our classifications. 

Table 1. Overview of the 15 genuine Herbig Ae/Be stars with 
their dominant spectral component (DC). 



Source 


DC 1 


Object 


DC T 


HD58647 


Si 


HD38120 


SP 


HD35929 


Si 


HD72106 


SP 


HD37258 


Si 


HD85567 


SP 


HD36917 


P 


HD95881 


SP 


HD37411 


P 


HD142527 


SP 


HD36112 


SP 


HD 190073 


SP 


HD37357 


SP 


HD244604 


SP 


HD37806 


SP 







1 Si : Silicate dominated, P : PAH dominated, SP : Silicates and PAHs 

To avoid the problem of separating the PAH bands from other 
dust components, we focus here on the analysis of the two PAH 
dominated spectra in the sample. In forthcoming papers the spec- 
tra dominated by silicate- or with mixed silicate-PAH emission 
features will be analysed (Bouwman et al. in prep.; Boersma et 
al. in prep.) 

The Herbig Ae/Be stars associated with these two spectra 
are HD36917 and HD3741 1. Both stars are located south of the 
Orion belt in the OB lc association, of w hich HD369 1 7 is a spec- 
troscopic binary dLevato & Abtl [l976). The shape of its SED, 
with its drop between 2-20 /im and rise beyond 20 /mi, sug- 
gest a cleared (inner) region in its circum-stellar disk (planet?). 
The star HD3741 1 is a visual binary, showing no emission lines. 
The optical spectrum of the binary indicates a late type (K-M) 
companion (Van den Ancker, unpublished). Table [2] summarises 
the astrometric data and Table[3]the available IRS data. FigureQ] 
presents the fully reduced Spitzer SL and SH spectra. 

Table 3. Available IRS data on HD36917 and HD3741 1. 







Module 








Object 


SL 1 


SH 1 


LH 1 


PU-blue 2 


# 3 


HD36917 
HD37411 


V(2) 
V(2) 


V(2) 
VO) 


V(4) 


V 


18 
22 



1 : Number of cycles in parentheses, 2 : Peak-up blue image (13.5 • 
18.7 /mi) 3 : Total number of spectra 



In order to derive the individual PAH profiles, a continuum 
is established. To facilitate comparison with previous studie s 
we adopt the procedures as outlined by iPeeters et al.l (|2002), 



2.0 r 



1.0 



0.8 - 



0.6 

0.5 
2.0 



t! 




0.25 
0.20 
S 0.15 
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6 8 10 12 14 
wavelength [micron] 


1 1 1 1 1 1 









HD36917 - 

-l 




/ 



HD37411 - 

9 10 11 12 13 14 15 20 
wavelength [micron] 

Fig.l. Fully-reduced SL and SH spectra for HD36917 and 
HD37411. The arrows indicate the positions of the PAH bands 
of interest. The inset displays the residual (SL) spectrum when 
approximating the continuum with a single spline. 



Ivan Diedenhoven et alj (|2004|) . and iHonv et all d200ll) . assum- 
ing PAHs are not significantly contributing to the continuum. In 
accordance with these procedures, the continuum for the 6.2 /im 
region is approximated by a single spline; for the 7-9 /im re- 
gion a general continuum and an additional plateau component 
are defined, where the general continuum is approximated by a 
single spline, with points selected between 5-6 /mi, 9-10 /im 
and near 7 /im. With additional points near 8.3 /im the plateau 
component is fixed. The continuum in the 11.2 /im region is ap- 
proximated by a single spline. Points between 9 - 10.5, 14.5 - 
15.5 and near 11.8 and 13.1 /mi are selected to define a single 
spline to approximate the continuum in the 12 /im region. 



3.1. Results 

Table [4] summarises the identified PAH bands and plateaus in 
the continuum subtracted spectra. The spectra of the two sources 
show all of the well-known IR emission features at 6.2, 7.7, 8.6, 
and 11.2 /im and, in addition, in some sources, the weak bands 
at 5.25, 5.7, 6.0, 12.0, 12.7, and 16.4 /im and the broad plateaus 
underneath the 6.2, 7.7, and 15 /im regions. Table |5]presents the 
spectral properties of the major PAH bands. These properties 
were determined fitting a single Gaussian to the band. 

Apart from PAH band emission also other, dust and non- 
dust components are identified. For inst ance, the aliphati c hy- 
drocarbon features at 6.85 and 7.25 /mi dSloan et alJl2005l and 
refs therein) in the SL spectra and the Nen, H2 S(l) and S III 
lines at, respectively, 12.8, 17.1, and 18.7 /im in the SH spectra 
of HD36917. Between 9.0 - 13 /im, a broad emission compo- 
nent, reminiscent of warm silicate dust, is visible, particularly in 
HD3741 1. In our analyses of the emission feature around 1 1 /im 
we subtracted this component using a spline fit (see Fig. 4). 



4 C. Boersma et al.: Chemical evolution in Herbig Ae stars 

Table 2. Astrometric data for HD36917 and HD3741 1. 



Object 


a [2000] 


<5 [2000] 


Region 


rf[pc] 


Ref. 


Sp. Type 


Ref. 


A v [mag] 


HD36917 


05h34m46.98s 


-05d34ml4.6s 


Orion OBlc 


510 


1 


B9.5V + A0V 


2 


0.49 


HD37411 


05h38ml4.51s 


-05d25ml3.3s 


Orion OBlc 


510 


1 


AO + K-M 


3 


0.00 



l: ldeZeeuw et"aT] ( ll999h ; 2: lLevato & AbtUl976l) : 3: lGrav & CorbalrvMl998h : 
Table 4. Detected PAH band components in HD36917 and HD3741 1. 



Source 


5.25 


5.7 


6.0 


6.2 


7.6/7.8 


8.6 


11.0 11.2 


12.0 


12.7 


13.5 


16.4 


Plateau 


HD36917 


V 


V 


V 


V 


V 


V 


V 


V 


V 


V 


V 


V 


HD37411 


V 




V 


V 


V 


V 


^ V 


9 


V 









Table 5. Spectral properties of the major PAH bands. 







HD36917 HD37411 


6.2 


position 1 
FWHM 1 
EW 2 
F 3 


6.22 6.25 
0.27 0.15 
0.028 0.049 
2.7 2.37 


7.7 


position 
FWHM 
EW 
F 


7.63 7.78 
0.42 0.43 
0.092 0.12 
6.3 3.63 


8.6 


position 
FWHM 
EW 
F 


8.58 8.64 
0.33 0.25 
0.022 0.036 
1.26 0.97 


11.2 


position 
FWHM 
EW 
F 


11.29 11.21 
0.31 0.38 
0.12 0.053 
3.74 0.94 


11.2 4 


position 
FWHM 
EW 
F 


11.32 11.22 
0.42 0.29 
0.11 0.063 
4.4 1.1 


12.7 


peak 
FWHM 
EW 
F 


12.72 12.71 
0.34 0.33 
0.026 0.013 
0.67 0.16 



1 : [fim], 2 : Equivalent width : J [(band - continuum) / continuumjdA 
3 : [10~ 15 • W • m" 2 ], 4 : Hi-res profile. 



variations in the 3.3 and 11.2 fim bands are studied. However, 
the observed variation are shown to be more modest and the ob- 
ject correlation not as tight than those observed for the 6-9 fim 
region. 

The observed diversity in peak position, profile, and rela- 
tive strength of the PAH bands are ascribed to global changes 
in the physical and chemical characteristics of the emitting 
PAH family. More specifically, variations in the peak posi- 
tions of the 6.2 fim band are thought to in dicate incorporation 
of nitrogen in the aromatic ring structure (Pee ters et al.l 120021 
lBauschlicherfl2"002l The 7.6/7.8 /zm bands are also likely to 
be related to chemical modifications, based on the strong cor- 
relation with the 6.2 fim variations. In contrast, variations in 
the ratio of the C-H modes to the C-C modes (e.g., 11.2/6.2 
fim) are attributed to va riations in the charge s tate of the 
emitt i ng PAH molecule s dAllamandola et all 1 19991: iHonv et al] 
120011; iBakes et alj|2001l) . Emission between 10.9 and 11.1 fim 
is attributed to the out-of-plane bending vibrations of solo- 
CH groups on the periphery of moderately sized (-100 car- 
bon atoms) PAHs dHudgins & Allamandoialll999l) . In Fig. |2] - 
[4] the derived profiles from the two Spitzer sources are com- 
pared to the characteris t ic pro files from lPeeters et al.l d2002l) and 
van Diedenhoven et al. (2004). 

Tabl e 6. The spe ctral characteristics for the classe s from 
IPeeters et all d2002l) and Ivan Diedenhoven et al.l d2004l) . X x in- 
dicates the peak positions, in fim. 



3.2. Profiles 

The apparent variations in peak position, width and wing shapes 
between the derived profiles is much larger than could be caused 
by the uncertainty in the adopted continuum. Studies of the de- 
tailed profiles of the PAH bands in the ISO SWS spectra of a 
large sample of stellar sources, planetary nebulae (PNe), reflec- 
tion nebulae, Hn regions and galaxies have revealed that the 
bands in the 6 - 9 fim range show strong variation in peak posi- 
tion and profile (IPeeters e t al. 2002). Those studies showed that 
the variations in the PAH bands correlate with object type. All 
ISM-like sources belong to group A characterised by a "6.2 //m" 
band peaking at 6.2 /jm and a 7.7 fim band peaking at 7.6 fim, 
and while the isolated Herbig Ae/Be stars, along with a few post- 
AGB stars and most PNe, belong to group B, which are charac- 
terised by a "6.2" fim band peaking around 6.28 fim and a 7.7 
fim band peaking at 7.9 fim. Finally, two post-AGB stars make 
up group C, characterised by a "6.2" fim band peaking at 6.3 fim 
and with no apparent 7.7 fim band, but a broad "8.22" fim fea- 
ture (Fig. [2]- El Table In Ivan Diedenhoven et all d2004l) . the 







The 5 - 9 fim region 






Class 


Object 


^6.2 


^7.7 1 


^8.6 


31 


ISM-like 


-6.22 


7.6/equal 


-8.6 


8 


non-ISM-like 


6.24 - 6.28 


"7.8" 


>8.62 


C 


non-ISM-like 


-6.3 


8.22 


none 






The 11.2 fim region 






Class 


Object 


^11.2 


FWHM, 


.2 


An.2 


ISM-like 


- 11.20- 11.24 


- 0.17 




A(B), 


2 uncorrelated 


- 11.20- 11.24 


- 0.21 




B] 1.2 


uncorrelated 


- 11.25 


-0.20 





1 : The "7.7" fim feature is classified by its dominant component, the 
7.6 and/or the 7.8 fim component. 



In Fig. |2l the 6.2 fim profiles are compared. The 6.2 fim 
profile of HD36917 resemblances that of class A best and for 
HD37411 that of class B. The profiles in both classes show a 
variation in the strength of the red tail, relative to the peak. In 
that respect the red tails of both sources fall off more rapidly 
than the chara cteristic profile that defines their class, but it is 
known to vary (IPeeters et al. 2002). 
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10.8 11.0 11.2 11.4 11.6 11.8 10.8 11.0 11.2 11.4 11.6 11.8 

wavelength [micron] wavelength [micron] 



Fig. 4. The derived 11.2 /mi profiles (black) after subtraction of the underlying continuu m and correction for the broad 10 /mi feature 
in HD37411. For comparison, the two characteristic profiles derived in the study by Ivan Diedenhoven et al. I (120041) are shown at 
matched resolution in grey. Left: The derived low-resolution (SL) profiles. Right: The derived high-resolution (SH) profiles. 



HD36917 




6.2 6.4 
wavelength [micron] 

Fig. 2. The derived 6.2 /mi profiles (black) after subtraction of 
the underlying continuum. For compar ison, the thre e chara cter- 
istic profiles derived in the study by iPeeters et alJ (120021) are 
shown at matched resolution in grey. 



8.0 8.5 9.0 9.5 

wavelength [micron] 

Fig. 3. The derived "7.7" /mi profiles (black) after subtraction 
of the underlying continuum. For comp arison, the th r ee ch arac- 
teristic profiles derived in the study by iPeeters et al.1 d2002l) are 
shown at matched resolution in grey. 



The variations in the "7.7" /mi feature (Fig. [3]) are particular 
interesting. Remarkably, HD36917 is best seen as a blend of a 
class A and class B profile. The "7.7" /mi band in the spectrum 
of HD3741 1 is a clear class B profile. 

The differences between the 1 1.2 /mi profiles are subtle and 
difficult to discern, despite the high-resolution profiles (Fig. 0). 
Again, HD36917 shows evidence for a blend of a class A and 
class B profile. The classification for HD37411 is hampered by 



the dominating continuum component between 9-13 /mi (inset 
Fig.[T|i. As described in Sect. 13. II this additional component has 
been approximated by a single spline and has been subtracted. 
The result is presented in Fig. |4] 

Both spectra show evidence for a 11.1 /mi feature. Many 
sources show this weak 's atellite' accompanying the 11.2 /mi 
feature dHonv et al.l 1200 lb . However, the 11.1 /mi feature in 
HD3741 1 is very strong relative to the 1 1.2 /mi band. This spec- 
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tral region is often confused by emission bands due to forsterite. 
In astronomical spectra, these often coincide with the 11.2 pm 
band. However, we note th at laboratory spectra show evidence 
for forsterite near 11.0 pm (Tam anai et alj |2006). The broad 10 
pm band, evident in Fig.Q] indicates the presence of warm sili- 
cate in the source and supports this suggestion. In any case, fo- 
cussing on the PAHs and guided by the high-resolution data, the 
profile is classified A. Table [7] summarises the assigned classes. 



Table 7. Assigned classes, see text for details. 
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Fig. 6. Correlation between the central wavelength of the 7-8 
pm band and the effective temperature of the exciting star. The 
5 Herbig stars used in th i s stud y are shown as black squares. 
Adopted from lSloan et all d2007l) 



4. Discussion 



ISloan et aD d2007h have identified systematic variations in the 
peak position of the 7.7 and 11.2 /mi band with the spec- 
tral type of the exciting star. We also recognise such varia- 
tions in the spectra of HD39917 and HD37411. We note that 
both s tars have similar spectr al type (~B9). Going back to ISO 
SWS ddeGraauw et all 1 19961) data dVan Kerckhovenll2002l) we 
have retrieved additional spectra of three well-studied and well- 
characterised Herbig Ae/Be stars, which provide additional in- 
sight in the possible factors driving these spectral variations. 

The first star is TY CrA, a dou ble-lined eclipsi ng binary 
consisting of a Herbig Ae/Be object dThe et alJll994l) . close to 
the zero age main sequence, and a late-type pre-main-sequence 
object. The object is located near one of the densest parts of 
the R Corona Australis star-forming re gion, embedded i n the 
bright reflection nebula NGC6726/6727 dHariu et alJll993l) . The 
infrared emission seen in the ISO spectrum does not origi- 
nate from TY CrA itself, but from the near by 'TY CrA bar' 
( Siebenmorgen et al. 20001: iGeers et al.ll2007l) . The second star 
is HD97048, it is located in the Chamaeleon I cloud of the 
Chamaeleon T association, which illuminates the bright reflec- 
tion nebula Ced 111. Spatial studies on the PAH emission done 
with the VLT's spectromete r and imager for the mid-infrared 
(VISIR; lLagageet all 120041) and ISO have revealed a circum- 
stellar (~ 0-5") and a more extended component (> 5.4"; 
Prusti et al.1 119941; ISiebenmorgen et alj 120001; Ivan Boekel et al.1 
2004 lHabart et alj|2006l) . Finally. HD100546 is a typical iso- 



lated Herbig Ae/Be star, which has no companions within 1500 
AU and is not associate d with the nearby dark cloud DC 292.6- 
7.9 dGradv etalj 12001). The mid-IR emission of this source 
clearly originates from the circumstellar protoplanetary disk, 
which has a gap of about ~10 AU with a 'puffed-up' inner rim 
dAcke & van den Anckerj|2006l) . 



4. 1 . The profiles 

Figure [5] compares the "7.7" //m feature profile^ of all five 
sources. The peak position of the 7.7 pm feature varies con- 

4 For details concernin g the ISO SWS spectra , such as data reduction 
and analyses we refer to Van Kerckhoven (2002). 



siderably in this sample from 7.6 pm all the way to 8.0 pm. 
Most striking is that the profile for HD97048 can be ob- 
tained by combining th e profiles from TY CrA and HD 100546 
dVan Kerckhovenl[2002h . The two Spitzer sources show similar 
behaviour. Indeed, the 7.7 /mi feature of HD36917 may be ob- 
tained by combining the class B profile of HD3741 1 with a class 
A-like profile (labeled pseudo) in Fig. [5] 

Table [7] shows that the systematic differences between the 
profiles of the three ISO sources are also analog to those be- 
tween the profiles of the two Spitzer sources. The spectral type 
of the five sources are similar (~B9), therefore, the variations 
in the derived profiles between these stars cannot be caused by 
diffe rences in the excit ing spectrum. 

ISloan et all d2007l) find a correlation between the central 
wavelength of the 7-8 pm feature and the effective tempera- 
ture of the exciting star. We have determined the peak position 
of the 7-8 pm band for the 2 Spitzer and 3 ISO sources in 
our study in the same manner as Sloa n et alTj 2007b and plotted 
our results with the data from ISloan et all "(120051 120071) in Fig. 
[6] Although there is no clear break in the c orrelation, a jump i s 
discernible around 10,000 Kelvin. Following Sloan et alj d2007l) . 
the diagram is interpreted to reveal a "7.7" pm band peaking at 
7.6 pm (class A) for T eff > 10,000 K, at 7.9 pm (class B) for 
10,000 > T eff > 7000 K and at 8.2 pm (class C) for T eff < 7000 
K. We note that this subdivision parallels a subdivision in source 
classification as with class A containing sources dominated by 
ISM material, class B comprising Herbig stars, and class C con- 
sisting of post-AGB objects. We point out that in this sample 
objects, like the PNe NGC 7027, with T eff = 200,000 K and be- 
longing to class B, are not shown. 

4.2. Extended PAH band emission 



IVan Kerckhovenl (|2002|) and iPeeters et all d2003l) note that the 
observed variations in the 7-9 pm region might well be due to 
differences in morphology of the sources, combining ISM ma- 
terial and/or disks. Most resounding in this respect is the source 
HD97048, which morphologically consists of a spatial combi- 
nation of circumste llar disk material and interstellar material 
dLagage et al.1 12004) and spectroscopically combines character- 
istics of classes A and B. 
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Fig. 5. Left: The normalised profiles of the "7.7" pm feature of HD97048, TY CrA, and HD100546, as obtained by the SWS 
instrument on board ISO. The dashed lines indicate the peak positions of the major PAH band components. In grey the combined 
profile for HD97048, constructed from TY CrA and HD100546, is presented. Note the quality of the match. Right: The normalised 
profiles of the "7.7" pm feature of HD36917 and HD3741 1, as obtained by the IRS instrument on board Spitzer. The profile labled 
pseudo is obtained by subtracting the profile of HD37411 from that of HD36917. The dashed lines indicate the peak positions of 
the major PAH band components. 



Linking back to the two Spitzer sources, this connects the 
detected PAH emission from HD3741 1 to the circumstellar disk 
and, particularly fascinating, the PAH emission from HD36917 
to, in part, the surrounding cloud and/or envelope (pseudo 
source) and, in part, the circumstellar disk. Intrigued by the 
notion of two different, spatially, separated PAH families in 
HD36917, we are motivated to investigate the spatial variation 
of the PAH emission coming from HD36917. 

The available 8 pm IRAC dFazio et all [2004) mosaic from 
part of the OBlc association is presented in Fig. [7] The image 
is centred on the position of HD36917 and the colour scaling is 
chosen to emphasise the emission from the envelope. The strong 
and rich morphology of the infrared background is clearly visi- 
ble. The IRAC data suggest diffuse extended nebulosity associ- 
ated with the source on scales up to at least ~ 1' (3 • 10 4 AU). 

Fig.[8]illustrates the use of the source profile extraction to de- 
compose both spatial components (see Sect.|2]l. In the left panel, 
part of the first order in the SL module of the spectral image 
of HD36917 is shown, demonstrating the extended nature of the 
1 1.2 and 12.7 ^m PAH bands in HD36917. The panel in the mid- 
dle displays the result (black) of simultaneous fitting a source 
profile and a 2 nd order polynomial (red) to the combined dither 



data (grey). A 2 order polynomial for the extended emission 
produces the optimal extraction, providing a reasonable approx- 
imation to the observed emission. A higher order polynomial re- 
sults in increased noise, while a lower order polynomial provides 
an inferior fit and underestimates the extended emission compo- 
nent. The half-width of the extended emission is about 6 pixels 
(~ 11" ; 6 • 10 3 AU). The panel to the right in Fig. [8] shows the 
cross-dispersion profile from 1 1.2 to 1 1.5 pm (black) compared 
with the calibration profile derived from standard stars (grey). In 
Fig. |9j the SL and SH spectra, corrected for local background 
emission, are presented. 

Remarkably, in both the SL and SH spectra the PAH sig- 
natures have largely disappeared. As expected, this procedure 
removes the forbidden atomic lines in the SH spectrum. In addi- 
tion, the 1 1.2, and 12.7 pm band and the plateau underneath 16.4 
pm have vanished. Obviously, the red end of the "7.7" pm fea- 
ture, the 11.2 and 12.7 pm PAH band emission originates from 
the surrounding cloud/envelope. It is clear that in this source 
there is no evidence for PAH emission associated with the di- 
rect (unresolved) circumstellar disk (~ 0.2" ; R ~ 50 AU). 

Note that, in contrast to the PAH features, the 16.4 pm band 
seems to be associated with the direct circumstellar environment 
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Fig. 7. 8 /^m IRAC mosaic from part of the OBlc association, 
centred on HD36917 and the colour scaling is chosen to empha- 
sise the emission coming from the envelope. Also shown is the 
IRS SL slit on the source. 



of the star. Although the band is located in a region where two 
orders overlap (12 and 13) and, therefore, the extended emission 
estimate is somewhat more uncert ain, forsterite does show some 
weak features near this position (Ko ike et"aTll2003l 16.3 //m). 
This actually suggests the presence of processed silicate mate- 
rial in the disk. Furthermore, spectral structure between 10.5 - 
10.9 fim appears in the background corrected spectra, most dis- 
tinctly seen in the SL s pectrum. This could be the weak PAH 
band first recognised by IWitteborn et"ai] ([l989) and later, e.g., 
iBeintema etafld 19961) . 

The absence o f PAH featur e s in th e spectra of the disk is 
not too surprising. Me eus et"ai] ([200 1) classified 46 spectra of 
Herbig Ae/Be stars, based on the shape of the SED in the mid- 
IR, in distinctive groups. Their group I sources all showed emis- 
sion due to PAHs, but their group II sources show no or only 
weak PAH emission. They link this difference to the geometry 
of the circumstellar disk; flaring vs. se lf-shadowed. The S ED of 
HD36917 makes it a Group II source (lHabart et al.ll2004 . con- 
sistent with this interpretation. 

The PAH emission in HD36917 originates from a region 
6 ■ 10 3 AU in size and, therefore, is likely associated with an 
extended envelope rather than a disk. Thus, although the spectra 
show evidence for a 'processed' PAH component peaking to- 
ward the source - thus associated with the source - the process- 
ing itself seems to occur on a scale-size, which is much larger 
than a star/disk system. Therefore, it has the appearance that the 
spectral variations can be attributed to the spatial structure of the 
region, although division into a clear disk and cloud component 
is not evident from our data. 

4.3. Origin of the variations 

The observed trend between the central wavelength of the 7.7 
- 7.9 jum feature and the spectral type of the irradiation star, as 




Fig. 8. Left: The normalised source profile of a section from the 
SL1 long-slit 2-D spectral image as observed toward HD36917. 
As a result of the angle of the order across the image, the pro- 
file centre moves in cross-dispersion. This results in a changing 
distribution of the profile signal across pixels, hence the appar- 
ent lower than unity profile at certain wavelengths. The pres- 
ence of a spatially extended component near the 11.2 and 12.7 
jum PAH bands in HD36917 is clearly visible. Middle : The si- 
multaneous fit (black) of a source profile and a 2 2n order poly- 
nomial to the combined dither data (grey). Right: The average 
cross-dispersion profile from 1 1.2 to 1 1.5 /urn is compared with 
the calibration profile derived from the standard stars. For the 
extraction the calibration cross-dispersion profile is fit in com- 
bination with a 2 nd order polynomial to decompose the spatial 
components. 



shown in Fig. [6] shows a jump near 10,000 Kelvin. The spec- 
tral variat ions in the "7.7" /xm feature reflect chemical modifica- 
tions and Sloa n et al. I (120071) suggested that these modifications 
are likely driven by UV processing. But hot class B sources, such 
as NGC 7027, seem to indicate that other source characteristics, 
besides effective temperature, may be of importance too. Such 
characteristics could include local density or history of the car- 
riers. Now, the jump near 10,000 Kelvin seems to indicate that 
for Herbig Ae stars these other factors are of consequence. In 
this study the spectral variations are linked to PAHs in spatial 
structure, promoting differences in UV processing. Note that the 
extended nature of HD36917 lends itself particularly well to a 
spatial-spectral study and may well settle the factors involved in 
the chemical processing of PAHs. 

4.4. Implications 

The chemical modifications driving the observed variations in 
Herbig Ae stars may well be due to a transition from sta- 
ble aromatic structures in the cloud (class A), to more la- 
bile aliph atic-like structure s in the circumstellar environment 
(class B) (Slo an et"aT]|2007l) . However, whereas post-AGB ob- 
jects (class C) and ISM sources (class A) lend themselves to 
an evolutionary interpretation - from labile aliphatic hydrocar- 
bons in the benign conditions of post-AGB objects to stable aro- 
matic hydrocarbons in the harsh environment of the ISM - it 
would be unlikely for the aromatic structures to regenerate their 
aliphatic side-groups when going from the ISM to the proto- 
planetary disks. Rather, we interpret the observed spectral varia- 
tions as reflecting the presence of an active chemical balance - in 
all sources - between hydrogenation, carbon reactions building 
(aliphatic) hydrocarbons, and UV photons breaking them down. 
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Fig. 9. The local background corrected SL (bottom) and SH (top) 
spectra for HD36917, as obtained with the decomposition. The 
corrected spectra are grey, the uncorrected spectra are dotted 
grey and the background estimates are black. The boxes show 
a few interesting spectral regions. 



5. Summary and conclusions 

From a sample of 15 Herbig Ae/Be stars we analysed the two 
PAH dominated Spitzer IRS spectra. The derived profi l es hav e 
been classified according the scheme of iPeeters et alJ ((2002). 
Comparison of these profiles with profiles obtained using ISO's 
SWS from three (spatially) well studied objects indicate the 
presence of two dissimilar, spatially separated, PAH families for 
HD36917. 

The analysis presented here shows again, based on the vari- 
ation in band profiles, that PAHs in subsequent, evolutionary 
linked stages of star formation , are different from those in the 
general ISM dHonv et al.ll200ll) . While previously the required 
chemical proces sing of the PAHs was placed in the protop lane- 
tary disk phase jVan Kerckhovenll2002l: IPeeters et aT1l2003l) . the 
evidence presented here, shows that such chemical changes may 
already occur in the (collapsing?) natal cloud. The supporting 
evidence for this is that the detected PAH emission cannot be 
associated with the (unresolved) disk and must thus be associ- 
ated with the circumstellar (natal) cloud. Furthermore, the pres- 
ence of the two distinct PAH families suggests active chemistry. 
Further studies on the spatial and spectral variations in the PAH 
spectra of these sources may be helpful to disentangle the factors 
that stimulate the processing of PAH molecules. 

Future observation, using (high-resolution) spectral mapping 
techniques, should provide a better handle on the spatial varia- 
tions of the PAH emission/profiles coming from and surrounding 



HD36917 and possibly link them to the changing radiation field 
as one moves out from the central star. Additional laboratory 
studies are required to explain the transition from class A PAH 
profiles in the circumstellar (natal) cloud to class B profiles in 
disks. 
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